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Предлагается новый алгоритм для формирования случайных фазовых экранов, используемых при численном решении задачи 
распространения лазерного излучения в турбулентной атмосфере. Моделирование экранов выполняется с учетом временных 
флуктуаций показателя преломления среды, что позволяет избавиться от периодичности, возникающей при поперечном смеще- 
нии экранов, имеющих конечные размеры. Таким образом, становится возможным формировать меняющиеся во времени турбу- 
лентные экраны неограниченной (бесконечной) протяженности. Новый подход находит свое применение в задачах моделирова- 
ния долговременных экспериментов в турбулентной атмосфере. Проведено сравнение численных и реальных экспериментов. 


Введение 

Исследования проблем, возникающих при рас- 
пространении лазерного излучения в турбулентной 
атмосфере, получили свое развитие еще в 70-х гг. 
прошлого века. К этому времени относятся теоре- 
тические результаты по выявлению основных эф- 
фектов и закономерностей в этой области. В на- 
стоящее время интерес исследователей привлечен 
к средствам компьютерного моделирования. 
Прежде всего, это связано с работой научных орга- 
низаций мирового сообщества по построению раз- 
личных лазерных систем, работающих в условиях 
реальной атмосферы. Построение каждой системы 
такого уровня связано с обязательными тестовыми 
испытаниями как системы в целом, так и ее состав- 
ных частей. Объем проводимых натурных экспери- 
ментов желательно минимизировать, т. к. постро- 
ение каждого варианта системы связано с больши- 
ми финансовыми затратами. Именно поэтому ме- 
тоды компьютерного моделирования становятся 
все более популярными. Благодаря сравнительно 
невысокой стоимости их реализации, возможности 
оперативного внесения изменений и отработки 
большого числа вариантов, развитию данных мето- 
дов уделяется все большее внимание. Ставятся за- 
дачи создания численных алгоритмов моделирова- 
ния, наиболее адекватно отражающих реальные ат- 
мосферные процессы, в том числе и флуктуации 
показателя преломления. 

При моделировании реальных лазерных систем 
важным требованием является выполнение расче- 
тов для достаточно продолжительного времени их 
работы. Это время может быть настолько велико, 
что возникает необходимость учета перемещения 
неоднородностей атмосферы ветром на расстоя- 
ния, существенно превышающие поперечные раз- 
меры пучка. 

Моделью турбулентных неоднородностей в на- 
иболее распространенных на сегодняшний день 
методах моделирования является один или нес- 
колько случайных турбулентных экранов [1—3]. 
Существующие методы формирования экранов 
ограниченной протяженности (порядка 2.. .20 ра- 
диусов пучка) непригодны для тех задач, где требу- 


ются экраны, существенно превышающие размеры 
радиуса пучка. В настоящее время авторам извест- 
на единственная работа, описывающая построение 
непериодических экранов произвольной протя- 
женности [4]. Особенностью данного подхода яв- 
ляется то, что он опирается на гипотезу Тейлора о 
замороженности турбулентности и требует боль- 
ших вычислительных затрат. 

В настоящей работе предлагается алгоритм 
формирования бесконечного турбулентного экра- 
на на основе идеи о его временной эволюции, что 
позволяет учесть дополнительную временную 
трансформацию распределения турбулентных нео- 
днородностей за счет флуктуаций скорости ветра. 

1. Алгоритм построения изменяющегося 

во времени экрана 

Традиционным подходом к решению задачи 
распространения лазерного излучения в турбу- 
лентной среде является метод расщепления по фи- 
зическим факторам [1]. При этом дистанция г (рас- 
стояние между источником и приемником излуче- 
ния) разбивается на N равных шагов Аі, а турбу- 
лентные неоднородности моделируются последо- 
вательностью случайных фазовых экранов. На 
каждом шаге влияние флуктуаций показателя пре- 
ломления атмосферы учитывается в приближении 
геометрической оптики. В этом случае функция 
корреляции флуктуаций фазы связана с трехмер- 
ным пространственным спектром флуктуаций по- 
казателя преломления соотношением [2]: 

В А г,Лг) = ||с/кФ х (к:,Дг) ехр(/кг) = 

= 2пк 2 Аг^ с/кФ п (к,Дг)ехр(жг), ( 1) 

где В 5 , - функция корреляция и спектр флукту- 

аций фазы, Ф„ - спектр флуктуаций показателя 
преломления, г ={х,у} - координата в плоскости, 
перпендикулярной оси распространиения, к -про- 
странственная частота в спектральном простран- 
стве, к - волновое число. При этом шаг Д^ должен 
быть, с одной стороны, достаточно мал, чтобы 
можно было использовать приближение геометри- 
ческой оптики, с другой стороны, достаточно ве- 
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лик, т. к. ур. (1) получено в приближении, что ди- 
станция много больше продольного простран- 
ственного масштаба корреляции флуктуаций пока- 
зателя преломления. 

Представляя случайный фазовый экран в виде 
преобразования Фурье (метод случайных спек- 
тральных амплитуд): 

5(г) = Лг/к:у4(к:)ехр(/к: г ) (2) 

и определяя, каким условиям должна удовлетво- 
рять спектральная амплитуда А(к), чтобы корреля- 
ционная функция для А(г) совпадала с Б/г,Дг), по- 
лучаем, что квадрат модуля спектральной амплиту- 
ды должен удовлетворять условию [3] 

|Л(к )| 2 =В$(к,Аі), 

где В ; (к) = — [ [ т/г/? ѵ (г)схр(-жг), т. е. аргумент 

случайной спектральной амплитуды является де- 
терминированным и его значение определяется ви- 
дом спектра флуктуаций фазы плоской волны, рас- 
пространяющейся в турбулентной атмосфере. Слу- 
чайность спектральной амплитуды определяется 
заданием ее случайной фазы. Выражение для слу- 
чайного фазового экрана имеет следующий вид: 

Д(г) = РТ {уІФ,(к) ехр[/ф(к)]} , (2) 

где РТ — оператор преобразования Фурье, ф(к) - 
случайный дельта-коррелированный фазовый 
множитель, равномерно распределенный в интер- 
вале [0, 2л]. 

Для решения динамических задач, связанных с 
поперечным движением турбулентных неоднород- 
ностей, обычно используется гипотеза заморожен- 
ное™ турбулентности 8(т,і)=5(г—Ѵ(,0), где У - по- 
перечная к оси распространения скорость ветра, 
1 - время. Используя свойство преобразования Фу- 
рье о сдвиге оригинала, из (2) получим следующее 
выражение для эволюции фазового экрана: 

У(г, г) = РТ {^Ф Дк) ехр[/ф(к)]ехр[ ік V?]}. (3) 

Отметим, что традиционная процедура форми- 
рования фазовых экранов на основе гипотезы Тейло- 
ра о замороженной турбулентности позволяет учесть 
перенос турбулентных неоднородностей ветром, но 
не позволяет учитывать флуктуации этих неоднород- 
ностей, присутствующие в реальной атмосфере. 

Сущность предлагаемого нами метода заключа- 
ется в том, что стандартный подход моделирования 
турбулентных экранов на основе гипотезы Тейлора 
(3) модифицируется к виду 

5(г , ( + ДО = РТ{^Ф 5 (к) /(< г + Д*,к) ехр[ ік ѴДг] }, (4) 

где изменение комплексной функции / во времени 
задается следующим рекуррентным выражением: 


1(1 + Аг, к) — р /(Ок) + V 1 -Р 2 ехр[кр(ок )], (5) 

где р=ехр(-А//г), А/ - шаг дискретизации задачи 
по времени, т - параметр, характеризующий ско- 
рость эволюции экрана (скорость флуктуаций по- 
казателя преломления среды), ф(1,к) - случайный 
дельта-коррелированный во времени (по 1) и в про- 
странстве спектральных амплитуд (по к) фазовый 
множитель, равномерно распределенный в интер- 
вале [0, 2 л]. Начальное выражение для функции / 
выбирается в виде 

/( 1 = 0,к) = ехр[іф(і = 0,к )]. (6) 

Выражения (4), (5) удовлетворяет необходимым 
условиям равенства нулю среднего значения и сох- 
ранения дисперсии случайной величины /[4]. Еще 
одним важным преимуществом такого подхода яв- 
ляется то, что он несущественно увеличивает об- 
щее время построения экрана, которое главным 
образом определяется временем вычисления бы- 
строго преобразования Фурье. 

Очевидно, что если не использовать рекуррент- 
ное выражение (5), а в (4) подставить начальное 
значение функции /в виде (6), то выражение (4) бу- 
дет совпадать с традиционным подходом, соответ- 
ствующим гипотезе Тейлора. 

Использование выражений (4), (5) позволяет 
наряду с переносом турбулентных неоднородно- 
стей средним ветром учесть дополнительную вре- 
менную трансформацию распределения турбулент- 
ных неоднородностей за счет мелкомасштабных 
флуктуаций скорости ветра. Под мелкомасштабны- 
ми флуктуациями подразумеваются пространствен- 
ные неоднородности распределения скорости ве- 
тра, размер которых меньше или порядка размера 
пучка. Коэффициент р=ехр(-Д?/г) задается в ин- 
тервале 0</кі и позволяет регулировать относи- 
тельный вклад случайных флуктуаций ветра на 
трансформацию неоднородностей атмосферы во 
времени по сравнению с их средним переносом. В 
частности случай р= 1 (т=оо) соответствует гипотезе 
Тейлора. При т=0 имеем экран, полностью обно- 
вляемый на каждом шаге по времени. 

Отметим, что Ѵ=0 соответствует отсутствию 
среднего ветра или совпадению его направления с 
направлением оси распространения пучка. В этом 
случае средний ветер не вносит вклада во времен- 
ные изменения неоднородностей, и данные изме- 
нения происходят исключительно за счет случай- 
ных флуктуаций скорости ветра с характерным 
временем, определяемым параметром г. 

Одним из важных достоинств предлагаемого 
подхода, связанного с использованием формул (4), 
(5), является то, что он позволяет получать неогра- 
ниченный и непериодический фазовый экран для 
долговременного моделирования движущейся тур- 
булентной среды. 
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2. Исследование применимости алгоритма 

к атмосферным задачам моделирования 

Для исследования применимости предлагаемо- 
го алгоритма формирования турбулентных экранов 
к атмосферным задачам было проведено сравнение 
результатов работы алгоритма с эксперименталь- 
ными данными. Эксперименты проводились со- 
трудниками ИОА СО РАН А. Л. Афанасьевым и 
А.П. Ростовым. Излучение распространялось на 
фиксированной во времени горизонтальной атмо- 
сферной трассе длиной 300 м на высоте 2,5 м от по- 
верхности земли. Длина волны излучения равня- 
лась 1,06 мкм, радиус пучка 12 см, радиус приемной 
апертуры 12 см. Среднее значение поперечной к 
трассе скорости атмосферного ветра изменялось от 
эксперимента к эксперименту в интервале от 0, 1 до 
2 м/с, среднеквадратичное отклонение скорости 
ветра от среднего значения о составляло от 0,2 до 


1,2 м/с. Время между последовательными наблю- 
дениями составляло 40 мс. Измерялись временные 
зависимости флуктуаций центра тяжести оптиче- 
ского изображения по осям хи у. Типичный вид 
данных зависимостей для различных среднеква- 
дратичных отклонений скорости ветра от среднего 
значения представлен на рис. 1. Видно, что в реаль- 
ном эксперименте присутствуют два характерных 
масштаба флуктуаций. 

Если построить автокорреляционные функции 
данных зависимостей 

у (А?) = — ^ — ((-ДО - < х>)(х{1 + Аі)- < х> )), 

то можно заметить, что степень корреляции в случае, 
соответствующем высоким флуктуациям скорости ве- 
тра а=0,5 м/с, падает быстрее, чем в случае, соответ- 
ствующем низким флуктуациям о=0,3 м/с (рис. 2). 



Рис. 1. Типичный десятисекундный интервал зависимостей флуктуаций х-координаты центра тяжести изображения для 
о=0,3 и 0,5 м/с. Среднее значение поперечной к трассе скорости атмосферного ветра составляло 0,44 м/с 



Рис. 2. Автокорреляционные функции флуктуаций х-координаты центра тяжести изображения, построенные по зависимостям, 
изображенным на рис. 1 
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Рис. 3. Моделирование замороженной турбулентности ( т=<уо) 



Рис. 4. Включение временной эволюции экрана т=1 с, а=0,3 м/с 


Нами проведен численный эксперимент для та- 
ких же параметров трассы и излучения для различ- 
ных значений параметра г. Результаты для т=со, 
что соответствует использованию традиционного 
алгоритма с ограниченным замороженным экра- 
ном, представлены на рис. 3. 

Как видно из сравнения рис. 1 и 3, использова- 
ние гипотезы Тейлора о замороженной турбулент- 
ности не позволяет наблюдать мелкомасштабные 
флуктуации центра тяжести изображения. Кроме 
того, использование ограниченных фазовых экра- 
нов приводит к тому, что в эволюции флуктуаций 
координат центра тяжести изображения наблюда- 
ется периодичность, в данном случае Т= 22 с. 

Как уже отмечалось ранее, предлагаемый в дан- 
ной работе алгоритм построения экрана позволяет 
избавиться от периодичности и учесть временную 
трансформацию распределения турбулентных нео- 
днородностей за счет мелкомасштабных флуктуа- 
ций скорости ветра. Положим г=2 с, что соответ- 
ствует снижению коэффициента корреляции экра- 
нов до уровня 0,2 за время 2,5 с. Время 2,5 с для за- 
данной скорости ветра и размеров экрана эквива- 
лентно времени полного проноса экрана вдоль оси 
х. Таким образом, для заданных начальных данных 
при г=2 с будет получен непериодический беско- 
нечный фазовый экран. Видно, что в отличие от 
случая т=со во временной зависимости флуктуаций 


центра тяжести изображения появляются мелко- 
масштабные флуктуации и одновременно пропада- 
ет периодичность (рис. 4). 

Отметим, что характер мелкомасштабных флу- 
ктуаций центра тяжести изображения, полученных 
в численном моделировании при т=1 с отличается 
от эксперимента, соответствующего о=0,5 м/с. Ис- 
следуя различные значения параметра т, можно по- 
добрать для него такое значение, при котором ре- 
зультаты численного эксперимента максимально 
приближены к результатам реального эксперимен- 
та. Анализ полученных результатов будем прово- 
дить на основе сравнения автокорреляционных 
функций, определяемых характером флуктуаций 
координат центра тяжести изображения. На рис. 5 
представлены автокорреляционные функции для 
результатов моделирования. 

Видно, что для наиболее адекватного модели- 
рования рассматриваемых натурных эксперимен- 
тов, соответствующих ст=0,3 и 0,5 м/с, подходят 
значения г=1 и 0,1 с, соответственно. Для модели- 
рования экспериментов с другими атмосферными 
условиями подбор соответствующего значения г 
также не представляет трудности. Кроме того, 
если необходимо моделировать изменение атмо- 
сферных условий в ходе эксперимента, можно 
варьировать значение г прямо в процессе модели- 
рования. 
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Рис. 5. Автокорреляционные функции для флуктуаций центра тяжести изображения при т =0,1; 1; 10 с. Результаты численного 
моделирования 


Заключение 

Для моделирования задачи распространения 
оптического излучения в атмосфере требуется учет 
временных изменений неоднородностей показате- 
ля преломления за счет мелкомасштабных флукту- 
аций скорости ветра. 

Предлагаемый в работе алгоритм построения 
фазовых экранов для моделирования влияния ат- 
мосферной турбулентности в предельном случае 
т=со совпадает с традиционным алгоритмом фор- 
мирования ограниченного экрана без учета вре- 
менной эволюции. Выбор параметра г позволяет 
варьировать скорость эволюции фазового экрана, 
учитывая тем самым флуктуации скорости ветра. 
Чем меньше значение г, тем быстрее эволюциони- 
рует экран и тем соответственно большему значе- 
нию флуктуаций скорости ветра оно соответствует. 

Если в ходе моделирования необходимо варьи- 
ровать атмосферные условия, что в рамках долго- 
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временных экспериментов представляется акту- 
альным, этого можно добиться соответствующим 
изменением параметра т. 

Еще одним немаловажным достоинством пред- 
лагаемого алгоритма формирования неограничен- 
ного экрана является то, что он не требует допол- 
нительных вычислительных ресурсов и его быстро- 
действие практически не уступает традиционному 
алгоритму формирования ограниченных фазовых 
экранов. 

Дальнейшее развитие предложенного в работе 
алгоритма к задачам атмосферной оптики авторы 
видят в определении количественных зависимо- 
стей используемого параметра г от параметров ат- 
мосферы и геометрии распространения излучения. 
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